Broadband diffraction management and self-collimation of white light in photonic 

lattices 
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We suggest a novel type of photonic structures where the strength of diffraction can be managed 
in a very broad frequency range. We introduce optimized arrays of curved waveguides where light 
beams experience wavelength-independent normal, anomalous, or zero diffraction. Our results sug- 
gest novel opportunities for efficient self-collimation, focusing, and reshaping of beams produced by 
white-light and super- continuum sources. We also show how to manipulate light patterns through 
multicolor Talbot effect, which is possible neither in free space nor in conventional photonic lattices. 
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It is known that periodic photonic structures can be 
employed to engineer and control the fundamental prop- 
erties of light propagation Q, Q ■ In particular, the beam 
refraction and diffraction can be modified dramatically, 
resulting in many unusual phenomena. For example, 
a beam can experience negative refraction in the direc- 
tion opposite to normal at the interface with a photonic 
crystal |2j, La- Lit LJ Additionally, the natural tendency 
of beams to broaden during propagation can be con- 
trolled through diffraction management [|J. Diffraction 
can be eliminated in periodic structures leading to self- 
collimation effect where the average beam width does not 
change over hundreds of free-space diffraction lengths • 
On the other hand, diffraction can be made negative al- 
lowing for focusing of diverging beams |8jl and imaging of 
objects with sub-wavelength resolution |3. Ho| . 

The physics of periodic photonic structures is governed 
by scattering of waves from modulations of the refractive 
index and their subsequent interference. This is a res- 
onant process, which is sensitive to both the frequency 
and propagation angle. Strong dependence of the beam 
refraction on the optical wavelength known as superprism 
effect was observed in photonic crystals Spatial 
beam diffraction also depends on the wavelength, and 
it was found in recent experiments 0, ^3 that the ef- 
fect of beam self-collimation is restricted to a spectral 
range of less than 10% of the central frequency. Such a 
strong dependence of the spatial beam dynamics on wave- 
length can be used for multiplexing and demultiplexing 
of signals in optical communication networks [Hi. flij . 
However, it remains an open question whether photonic 
structures can be used to perform spatial steering and 
shaping of beams emitted by white-light sources, such as 
light with supercontinuum frequency spectrum generated 
in photonic-crystal fibers and fiber tapers fl5t [uj . 

In this Letter, we suggest a novel type of periodic 
photonic structures designed for wavelength- independent 
diffraction management in a very broad frequency range, 
covering a spectral range up to 50% of the central fre- 
quency. We introduce the optimized periodic struc- 



tures where multicolor beams experience constant nor- 
mal, anomalous, or zero diffraction. This opens up 
novel opportunities for efficient self-collimation, focus- 
ing, and shaping of white-light beams. For example, in 
such optimized structures it becomes possible to manipu- 
late white-light patterns through multicolor Talbot effect, 
which otherwise is not feasible in free space or in conven- 
tional photonic lattices. 

We study propagation of beams emitted by a contin- 
uous white-light source in a periodic array of coupled 
optical waveguides [see Fig. Qfa)], where the waveguide 
axes are also periodically curved in the propagation di- 
rection [see examples in Figs.EJa) andOfa)]. In the linear 
regime, the overall beam dynamics is defined by indepen- 
dent evolution of complex beam envelopes E(x, z; A) at 
individual frequency components governed by the nor- 
malized paraxial equations, 
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—u[x-x (z)]E = 0, (1) 



where x and z are the transverse and propagation coor- 
dinates normalized to the characteristic values x s = Ifim 
and z s = 1mm, respectively, A is the vacuum wavelength, 
c is the speed of light, Uq is the average refractive index of 
the medium, v(x) = v(x + d) is the refractive index mod- 
ulated with the period d in the transverse direction, and 
xq(z) = xq(z + L) defines the longitudinal bending pro- 
file of the waveguide axis with the period L 3> d. When 
the tilt of beams and waveguides at the input facet is 
less than the Bragg angle at each wavelength, the beam 
propagation is primarily characterized by coupling be- 
tween the fundamental modes of the waveguides, and can 
be described by the tight-binding equations taking into 
account the periodic waveguide bending [l3.ll7|. 



dz 



^ r. 



■ ujx (z)n^ n , (2) 



where ^ n {z;uo) are the mode amplitudes, n is the waveg- 
uide number, w — 2imod/\ is the dimensionless fre- 
quency, and the dots stand for the derivatives. Coefficient 
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FIG. 1: Discrete diffraction in (a) straight waveguide array 
with period d = 9fim. (b) Coupling coefficient normalized 
to the coupling at the central frequency Co- (c-e) Evolution 
of beam intensity and output intensity profiles after 80mm 
propagation of a 3/im wide input beam for (c) A r = 580nm, 
(d) Ao = 532nm, and (e) \t ~ 490nm, which correspond to 
the points 'c', 'd', and 'e' in (b). Waveguide width is 3/xm 
and substrate refractive index is no — 2.35. 



C(uj) defines a coupling strength between the neighboring 
waveguides, and it characterizes diffraction in a straight 
waveguide array with xq = 0,0]. The coupling coeffi- 
cient decreases at higher frequencies [2(J and accordingly 
the beam broadening is substantially weaker at shorter 
wavelengths, see Figs.^b-e). 

We consider symmetric profiles of the waveguide bend- 
ing such that Xq(z) — f(z — z a ) for a given coordinate 
shift z a , where function f(z) is symmetric, f(z) = f(—z). 
Then, after a full bending period (z z + L) the beam 
diffraction is the same as in a straight waveguide array 
with the effective coupling coefficient 113 



C eff H = C(u)L- 



cos [wio(C)] ^C- 



(3) 



According to Eq. J2J, diffraction of multicolor beams is 
defined by an interplay of bending-induced dispersion 
and frequency dependence of the coupling coefficient in 
a straight waveguide array. We suggest that spatial evo- 
lution of all frequency components can be synchronized 
allowing for shaping and steering of multi-color beams, 
when effective coupling remains constant around the cen- 
tral frequency u)q, 



dC e s{u>) / du>\ 



0, 



(4) 



and we demonstrate below that this condition can be 
satisfied by introducing special bending profiles. 



First, we demonstrate the possibility for self- 
collimation of white-light beams, where all the wavelength 
components remain localized despite a nontrivial evolu- 
tion in the photonic structure. Self-collimation regime 
is realized when the diffraction is suppressed and the ef- 
fective coupling coefficient vanishes, C c g = 0. This ef- 
fect was previously observed for monochromatic beams 
in arrays with zigzag or sinusoidal [l^ bending pro- 
files, however in such structures the condition of zero 
coupling cannot be satisfied simultaneously with Eq. J2J, 
resulting in strong beam diffraction under frequency de- 
tuning by several percent jl^| . We find that broad- 
band diffraction management becomes possible in hy- 
brid structures with a periodic bending profile that con- 
sists of alternating segments [see example in Fig. Ela)], 
xq(z) — Ai{cos[2n z/z ] — 1} for < z < z , x (z) — 
A 2 {cos [2tt(z - z )/{L/2 - z a )} - 1} for z Q < z < L/2, 
and xq(z) = —xq(z — L/2) for L/2 < z < L. Effective 
coupling in the hybrid structure can be calculated analyt- 
ically, CeffM = C(^)2L- 1 [z J (Ci) + (L/2 ~ z )J {&)], 
where J m is the Bessel function of the first kind of the 
order m, £i = 2irAiui/ zq, and £2 = 27r^42w/ {L/2 — zq). 

We select a class of symmetric profiles of the waveg- 
uide bending to avoid asymmetric beam distortion due to 
higher-order effects such as third-order diffraction. Ad- 
ditionally, the waveguides are not tilted at the input, i.e. 
Xo(z = 0) = 0, in order to suppress excitation of higher- 
order photonic bands by incident beams inclined by less 
than the Bra gg ang le. The effect of Zener tunneling to 
higher bands [2l], |23 and associated scattering losses can 
be suppressed irrespective of the waveguide tilt inside the 
photonic structure by selecting sufficiently slow modula- 
tion to minimize the curvature xq (z) and thereby achieve 
adiabatic beam shaping. 

In order to realize broadband self-collimation, we 
choose the structure parameters such that £i(wo) = 
£1 ~ 2.40 and £2(^0) = £2 — 5.52 are the first and 
the second roots of equation Jo(0 = 0. Then, the 
self-collimation condition is exactly fulfilled at the cen- 
tral frequency ujq, C e ff(o;o) = 0, and simultaneously the 
condition of frequency-independent coupling in Eq. Q 
is satisfied for the following modulation parameters, 
A x = Ki6Ji(6)/2^(6Ji(6) -£Ji(£)>o]L/2, A 2 = 
-[Ji(t;i) / Ji(£,2)]Ai, and z = 27rw ^i/£i- As a result, 
we obtain an extremely flat coupling curve shown in 
Fig. EJb) where the point 'd' corresponds to the cen- 
tral frequency. In this hybrid structure not only the 
first derivative vanishes according to Eq. JIJ , but the sec- 



ond derivative vanishes as well, 



< 10- 



As a result, 



fiJa(Ci)Ji(6)-6J2(6)Ji(ei) 
the effective coupling remains close to zero in a very 
broad spectral region of up to 50% of the central fre- 
quency. We note that the modulation period L is a free 
parameter, and it can always be chosen sufficiently large 
to avoid scattering losses due to waveguide bending since 
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FIG. 3: Wavelength-independent diffraction in an optimized 
periodically curved waveguide array, (a) Waveguide bending 
profile with the period L = 40mm and (b) corresponding ef- 
fective coupling normalized to the coupling in the straight 
array at the central frequency Co = G (wo)- (c-e) Evolution 
of beam intensity and output intensity profiles after propaga- 
tion of two full periods for the wavelengths (c) A r = 580nm, 
(d) Ao = 532nm, and (e) A;, = 490nm, which correspond to 
points 'c', 'd', and 'e' in plot (b). 



FIG. 2: (a-e) Broadband self-collimation in an optimized 
waveguide array: (a) Waveguide bending profile with the pe- 
riod L = 60mm and modulation parameters A\ = 27/im, 
Ai = 42^im, zo — 18mm. (b) Effective coupling normal- 
ized to the coupling in the straight array at the central fre- 
quency Co = C(a>o). (c-e) Evolution of the beam intensity 
and output intensity profiles for different wavelengths marked 
(c) A r = 560nm, (d) Ao = 532nm, and (e) A(, = 400nm corre- 
sponding to marked points in (b). (f-h) Frequency-sensitive 
diffraction in array with the sinusoidal bending profile at the 
wavelengths corresponding to plots (c-e). 



the maximum waveguide curvature is inversely propor- 
tional to the period, max|xo(z)| ~ L^ 1 . Although the 
beam evolution inside the array does depend on the wave- 
length, the incident beam profile is exactly restored after 
a full modulation period, see examples in Figs. EJVe). 
Self-collimation is preserved even at the red spectral edge, 
where coupling length is the shortest and discrete diffrac- 
tion in the straight array is the strongest [cf. Fig. P^c) 
and Fig. ^c)] . The hybrid structure provides a dramatic 
improvement in the bandwidth for self-collimation effect 
compared to the array with a simple sinusoidal modu- 
lation, where beams exhibit diffraction under small fre- 
quency detuning, see Figs. |5f f-h). 

We now analyze the conditions for frequency- 



independent normal or anomalous diffraction that may 
find applications for reshaping of multicolor beams. In 
order to reduce the device dimensions, it is desirable 
to increase the absolute value of the effective coupling 
and simultaneously satisfy Eq. J3J to achieve broadband 
diffraction management. We find that Eq. J3J can be sat- 
isfied even in the simplest two-segment hybrid structure 
with zo = L/2 and A\ = (£/2iruio) L/2. Here a set of 
possible parameter values £ is determined form the re- 
lation Jo(£)/Ji(£) = Co^/Ciujo, where Co = C(u>o) and 
C\ = dC(u>)/duj\ ui=ujg characterize dispersion of coupling 
in a straight array. It is possible to obtain both normal 
and anomalous diffraction regimes for normally incident 
beams, corresponding to positive and negative effective 
couplings Ceff(wo) = CoJo(C) depending on the chosen 
value of £. For example, for the waveguide array shown 
in Fig. ^ at the central frequency u>q = 250 [correspond- 
ing wavelength is Ao = 532nm] coupling parameters 
are C ~ 0.13mm- 1 and C x ~ -0.0021mm- 1 . Then, 
constant positive coupling around the central frequency 
Ceff(wo) — 0.25Co is realized for £ ~ 6.47 and constant 
negative coupling C c s{^o) — — 0.25Co for £ ~ 2.97. 

We perform a comprehensive analytical and numeri- 
cal analysis, and find that a hybrid structure with bend- 
ing profile consisting of one straight (i.e A\ = 0) and 
one sinusoidal segment can provide considerably im- 
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FIG. 4: (a) Monochromatic Talbot effect in the straight 
waveguide array shown in Fig.^a): periodic intensity revivals 
every — 16.5mm of propagation for the input pattern 
{1, 0, 0, 1, 0, 0, . . .} and the wavelength Ao = 532nm. (b) Dis- 
appearance of the Talbot carpet in the straight array when 
input consists of three components with different wavelengths 
A r = 580nm, Ao = 532nm, and Xt = 490nm and equal inten- 
sities, (c) Multicolor Talbot effect in the optimized structure 
with wavelength-independent diffraction [see Fig. Half of 
the bending period L/2 = = 53.2mm is equal to the Tal- 
bot distance for the corresponding effective coupling length. 



proved performance if ujqCi/Cq > £, C rJi{ticr)/Ja{£,cr), 
where value £ cr ~ 5.84 is found from the equa- 
tion [Jl(£cr) + £,cr [Jo(Ccr) " J 2 (£cr)] /2] ]J (£cr) - 1] + 

£cr</i(£cr) = 0. Under such conditions, larger val- 
ues of positive effective coupling can be obtained 
in a hybrid structure with A\ = 0, Ai = 
[C 1 C off ( Wo )/2^C 2 J 1 (6)]£/2, z = [C cS (lo )/C q ]L/2. In 
this structure, the effective coupling at central frequency 

is OrffM = &cgJ 1 (6)/[e 2 CbJ 1 (6) +W0CJ. 

Example of a hybrid structure which provides strong 
wavelength- independent diffraction is shown in Fig.|3Ja), 
and the corresponding effective coupling is plotted in 
Fig. 03T>). The output diffraction profiles in this opti- 
mized structure are very similar in a broad spectral re- 
gion, see examples for three wavelengths in Figs. I^c-e). 
We note that the outputs at these wavelengths are sub- 
stantially different after the same propagation length in 
the straight waveguide array, as shown in Figs.^c-e). 

As one of the applications of the broadband diffrac- 
tion management we consider a multicolor Talbot effect 
which allows to manipulate white- light patterns. The 
Talbot effect, when any periodical monochromatic light 
pattern reappears upon propagation at certain equally 
spaced distances, has been known since the famous dis- 
covery in 1836 [23j. It was recently shown that the 
Talbot effect is also possible in discrete systems for 
certain periodic input patterns 20]. For example, for 
the monochromatic periodic input pattern of the form 
{1, 0, 0, 1, 0, 0, . . .}, Talbot revivals take place at the dis- 
tance = (2tt/3) [1/C(u]\, see Fig. Ufa) . 

Period of the discrete Talbot effect in the waveguide 
array is inversely proportional to the coupling coeffi- 
cient C{u>), which strongly depends on frequency, see 
Fig-Eb). Therefore, for each specific frequency Talbot 
recurrences occur at different distances 20], and peri- 



odic intensity revivals disappear for the multicolor input, 
see Fig. Etb). Multicolor Talbot effect is also not pos- 
sible in free space where revival period is proportional 
to frequency. Most remarkably, multicolor Talbot ef- 
fect can be observed in optimized waveguide arrays with 
wavelength-independent diffraction, see Fig.^c). In this 
example, we use the shape of structure with constant pos- 
itive diffraction shown in Fig. and choose half of the 
bending period to be equal to the period of the Talbot 

recurrences for the corresponding effective coupling in 
(2) 

this structure, L y T ' = (27r/3) [1/C e g(u))]. We note that 
the length of the straight segment is equal to the Tal- 
bot distance at central wavelength in the straight array, 
zq = , which explains partial revivals at the end of 
the straight segments visible in Fig. Uc). 

In conclusion, we have introduced a novel class of pho- 
tonic structures where diffraction can be engineered in 
a very broad frequency range. We have analyzed the 
optimized array of periodically curved waveguides where 
light beams experience wavelength-independent normal, 
anomalous, or zero diffraction, and predicted the multi- 
color Talbot effect which is not possible in free space and 
conventional waveguide arrays. Our results suggest novel 
opportunities for efficient self-collimation, focusing, and 
reshaping of beams produced by white-light and super- 
continuum sources. 
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